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1st AIAA Geometry and Mesh Generation Workshop 
3 AIAA High-Lift Prediction Workshop 
Fuselage 


Inboard and 


- Introduction Outboard Flaps 


Wing 
- Structured overset meshing methods and best practices 
using Chimera Grid Tools (CGT): AIAA Paper 2017-0362 


- Lessons learned 
- Meshing a family of grid systems at different resolutions 
- Grid quality checks 


- Summary and conclusions 
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STRUCTURED OVERSET MESHING USING CGT: 


METHODS AND BEST PRACTICES 
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A HIGH-LIFT CRM GEOMETRIC CONFIGURATIONS cS] 
Ames Research Center 


Full Flap Gap (coarse, medium, fine, extra fine) 


inboard flap— 


outboard flap— 


Partially-Sealed 
Flap Gap 
(medium only) 
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A OVERSET STRUCTURED GRID GENERATION So 
PROCESS AND SCRIPTING FRAMEWORK 
Main steps 
- Geometry processing 


- Surface grid generation: featured-based domain decomposition, 
grid point distribution, mesh fill 

- Volume grid generation: hyperbolic near-body, Cartesian off-body 

- Domain connectivity: grid points blanking, donor stencil search 

- Input parameters preparation for flow solver: 


boundary conditions, grid indices for component aerodynamic loads 


Develop script that reproduces entire process 


- Chimera Grid Tools Script Library (Tcl-based, 200+ macros) 


- Component scripts (fuselage, slat, flaps, wing) 
- Master script 


- Parameterized inputs 
- max stretching ratio (surface and volume) 


- surface grid spacing (max interior, at surface features) 
- volume grid wall normal spacing 


- min number of points on smallest feature 
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JA GEOMETRY PROCESSING S 


- Geometry definition files supplied: native CAD, STEP, IGES 
- Create starting point for grid generation script development 
- Unstructured surface triangulation (CART3D format) 
- Generated using ANSA software 
- Sufficient resolution at high curvature regions 
- Surface curves (PLOTSD format) 
- Generated using Chimera Grid Tools from surface triangulation 
- CAD edges including all surface features 
- Identify configuration characteristic lengths 
- component length scale 
- smallest feature size 
- gap size between components 


chord A 


gap size 


oN cusp / t.e. 
Pa thickness 


AA crip POINT DISTRIBUTION MESHING GUIDELINE 
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Mostly prescribed by High-Lift Prediction Workshop document 


Root 
spacing 


Nose 
spacing 


Leading edge 
— spacing 


Tail 
spacing Maximum interior Tip 
spacing and spacing 
stretching ratio 
Wall normal 


Streamwise spacing 
Number of points across at t.e. (Spacing 
blunt trailing edge/cusp-” continuity for 
J / O-topolog 
hg Trailing edge 
spacing 9——~ 


Wall normal 
stretching ratio 
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WORKSHOP PRESCRIBED MESHING PARAMETERS 
Reference spacing AS,,; = 3% mean aerodynamic chord 


Resolution Level 


# Points on trailing edge 


Span spacing at flap gap 
cap grids (x 10-2) * 


Max surface spacing 


* Not prescribed by workshop 


GRID QUALITY CHECKS 


All volume meshes are automatically checked for 


- Positive Jacobians as computed by target flow solver (OVERFLOW) 
- Self intersections with surface mesh 
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A INITIAL CURVES AND SURFACE GRIDS cS] 
Fuselage Features: Cockpit Window, Fairing, Wing 
Intersection 


-_ Cockpit Window 
Initial Curves (28) 


Geometry trimmed 
by wing 
intersection curve 


Fuselage Fairing 


Wing/Fuselage Intersection 


Fuselage Side of Wing/ 
Fuselage Collar Grid 


Surface Grids (13) 


L.. enter INITIAL CURVES AND SURFACE GRIDS 
Slat and Flap Features: L.E., T.E., Cusp, Root, Tip 


Slat Cap Grids 
(5 at each root 
and tip) 


YY 
YY 
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Slat Initial 


Curves (20) Flap Cap Grids 


(4 at each root 
and tip) 


Flap Initial 
Curves 
(2x12) 


Flap Surface 
Grids (2x10) 
Flap Gap 

Region 
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INITIAL CURVES AND SURFACE GRIDS 
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Wing Initial 
Curves (65) 
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L..... FLAP PARTIAL SEAL SURFACE GRIDS ] 


Re-use grids from full flap gap case for fuselage, slat, wing, and flaps 


Partial flap seal against fuselage Partial flap seal between inboard 
and outboard flaps 


Back and side wall cap split into two grids to avoid double 
concave corner => easier for hyperbolic volume mesh generation 
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SLICES OF FUSELAGE, SLAT, FLAP 
Ames Research Center VOLUME GRIDS 


Fuselage 


Slat Root 


Ebi, 
RON 
RAS 
WAAR 
RRS 
\' 


ie 
x 
SY 


S 
RS SS WX S 


LTT) 
CTT 


“ 
= 
o 
©) 
Wu 
= 
> 
mal 
O 
> 
© 
< 
= 
LL 
O 
Tp 
Wu 
o 
al 
a 


Wing/Fuselag 


Collar Gr 


® 
> 
fe) 
O 
Q 
= 
rT 


Side Wall 


Af 


ve 


—S 


A WALL-NORMAL GRID POINT DISTRIBUTION 


Fuselage and Slat 
1. Uniform spacing first two cells (As,,,,) 
2. Stretched region to outer boundary 


Stretched regions 
Shear layer region 


Wing and Flaps 
Need to resolve shear layer from preceding 
component for accurate drag prediction 
1. Uniform spacing first two cells (As,,,,)) 
2. Stretched region 
3. Shear layer region 
Uniform spacing = 100 x As,,,,, 
Thickness = 3 in. 
Distance from wall = 1.5 in. 
4. Stretched region to outer boundary 


Wing = 
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A OFF-BODY STRETCHED CARTESIAN S 
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- Cartesian box grid with uniform core and stretched outer layers 
- One box grid around fuselage volume grids 
- Three staggered box grids around slat, wing, flaps 


- Large stretched box grid encloses all 
smaller box grids and goes to far field 
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eee ee DOMAIN CONNECTIVITY os) 
Comparison of Two Approaches 


Chimera Components Connectivity Program (C3P) 


- Inputs: boundary conditions for each mesh, and component ID 
for each solid wall (low manual effort needed) 


- External process performed prior to running OVERFLOW flow solver 


OVERFLOW-DCF (DCF) 

- Inputs: boundary conditions for each mesh, X-ray map for each hole 
cutter, list of grids to be cut by each X-ray, constant offset distance 
for each hole cut instruction (significant manual effort needed) 


- Built into the OVERFLOW flow solver 
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Z > CONSTANT SPAN CUT THROUGH VOLUME MESH @D 
(Wing Trailing Edge, Flap Leading Edge Region) 
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(spatially variable offset) (constant offset) 
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VARIOUS VOLUME SLICES FROM C3P CONNECTIVITY 


Inboard Flap 


Fuselage 


M4 
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seeS 
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Constant-x cut at wing/ 
fuselage junction 


Constant-x cut across flap gap 
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Agri SCRIPT DEVELOPMENT TIME (MEDIUM MESH gy 


Two different domain connectivity methods/software 


Surface grid generation 


Volume grid generation | 450 | 66 | 


Input prep. (flow solver b.c., post-processing) | 34 | 45 | 
foal CST | 100 


Time (hr.) | % of Total 
Geometry processing 
Surface grid generation 


Volume grid generation | 450 | 56 | 


Input prep. (flow solver b.c., post-processing) | 34 | 389 
Total 80.2 | 100.0 
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GRID SCRIPT DEVELOPMENT FOR DIFFERENT LEVELS OF 
MESH RESOLUTION AND PARTIALLY-SEAL FLAP GAP 


Full flap gap coarse, fine, and extra-fine level grid systems 
- Created independently from the medium level system 
- Not a redistributed version of medium mesh 


Partially-sealed flap gap medium system 
- Created by copying fuselage, slat and wing grids, and some flap 
grids from full gap grid system 
- Only need to create grids for partial seals 


FlapGapGeometry |  sFullGap_ =| Partial Seal 
Resolution Level 


Grid script 17.75 * 
development time (hr.) 


* Additional development time beyond medium mesh script 


LZ o- 
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- Entire process performed on Linux Xeon desktop workstation 
- All timings include i/o 


[Fiap Gap Geometry (| —~S=SC~SulNGap——SSS—=*d én Seal 
Medium 
#GridsSSC*dSC 7 
50 
# Volume grid points (x 10°) | 24.1 63 
[Grid script devel. time (hr) | 100 | 686 |1775| 125 | 120 
[Grid seript exec. time (min) | 3.25 sas | 1.65 


73 
Connectivity (C3P) exec. wall 1.14 2.85 7.25 28.23 3.1 
time (min.), mem use (GB) (6) (13) (31) (81) (13) 
24 OpenMP threads 
Connectivity (DCF) exec. wall n/a n/a n/a 
time (min.) 
24 MPI ranks 


n/a = not attempted due to extra manual time needed to create special X-ray cutters 


| Fine 
| 1.02 
| 189.3 


Medium 
| 0.51 
| 65.4 
| 5.35 
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A OVERFLOW PRELIMINARY RESULTS S&S 
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(Tom Pulliam) 
Grid Sizes saab weirs re nas 
Coarse: 2.4M ee 
Medium: 6.5M 
Fine: 189M 
X-Fine: 565M 
Connectivity: C3P and DCF 


\ Cp On Surface 
-6.00 


Parameters: 
Mach = 0.2, Alpha = 8 deg, 
Re = 3.26M based on MAC 
3° order Roe, SA-RC-QCR2000 


HiLIft-CRM Workshop Grids 
SA-RC-QCR 2000, 3rd Order Roe, M = 0.2, a= 8.0° 
1.78 T | T T T 


C L @@ C3P | 
rs @ DCF 

L @ = Partial Gap Seal 

Ltt a 


OS 1.76 =] 
L- a a 


L 
0 Se-06 


ee oansity Contours” 


HiLift-CRM Workshop Grids 
C3P, SA-RC-QCR2000, 3” Roe, M = 0.2 0=8° 


TIN ©) N = # Grid pts 
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LESSONS LEARNED FROM MESHING EXERCISE 
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GENERATION OF FAMILY OF GRID SYSTEMS AT 
DIFFERENT MESH RESOLUTIONS 
(COARSE, MEDIUM, FINE, EXTRA-FINE) 


- Grid system at each resolution level is generated independently of 
each other starting from geometry definition 


- Different meshing parameters prescribed for each level (e.g., max 
stretching ratio, max interior surface grid spacing, grid spacing at 
surface features, number of points on t.e., volume mesh wall 
normal spacing) 


- Current practices do not provide automatic adjustments of 
marching distances and smoothing parameters 


- Significant grid script execution time at extra-fine level (> 0.5 hr) 


26 
JA PARAMETER ADJUSTMENTS AT cS] 
DIFFERENT LEVELS OF GRID RESOLUTIONS (1) 


Hard coded grid indices for medium mesh script 
Splitting locations defined by 
- Grid indices => faster one level (medium) mesh development 
- Grid coordinates or distance to reference point => 
slightly slower one level development but works for all levels 


Z 


Y 
Split locations Ley GO4 
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Ames Research Center PARAMETER ADJUSTMENTS AT 
DIFFERENT LEVELS OF GRID RESOLUTIONS (2 


Hyperbolic grid marching 
distances chosen to 
provide optimal overlap at 
coarse level 

(e.g., 5-point overlap for 5- 
point flow solver stencil) 
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- Too much overlap at fine ae 
and extra fine levels poe 
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“Ames Research Center PARAMETER ADJUSTMENTS AT 
DIFFERENT LEVELS OF GRID RESOLUTIONS (3) 


Finer grid spacing in concave corners in finer levels 
- Need to adjust smoothing parameters for hyperbolic marching 
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GRID QUALITY CHECK UTILITIES © - 
NOT CURRENTLY IN CGT 


Need min/max and distribution of grid attribute 
statistics => Histogram and color map display 


1. Distance of surface grid points to geometry 
definition (Native CAD, STEP, IGES) 


2. Distance to wall of first grid 
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3. Cell orthogonality 
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Jacobians and Cell Volumes 
Must pass 
1. Jacobian computed using same subroutine as in target flow solver 
OVERFLOW (all > 0) 
2. Self-intersection of volume grid points against surface grid (none) 


Mostly pass 
3. Cell volume using hexahedral decomposition into 6 tetrahedrons 
4. Stretching ratio (<= 1.2) 


Cell Volume Check 
Pd i Oe Sum 6 tets One or more tets 
Negative sum Negative vol 


Entity #Cells (J,K,L) Min Vol #Cells Min Tet Vol 


1 0 pa 
2 0 
3 0 2.526e04 9 -2.719e-06 
4 0 
- Cut into 2 prisms 5 0 9.463e-04 186 -2.195e-04 
P : 6 6 60 
- Cut each prism into 3 tets 7 0 4.672e-04 59 -2.892¢-04 
: 8 0 3.594e05 12 -4.785e07 
- Bad cell if 9 0 5.191e04 6 -1.232e07 
1. any tet volume <0 10 1 (26,38,28) -1.567e-03 216 -1.884e-03 
1 
2. sum of 6 tet volumes <0 12. 0 


OVERGRID Diagnostic 
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Orphan Points 


ivity 
Number, location and spread (OVERGRID) 


GRID QUALITY CHECK UTILITIES IN CGT (2) 
Domain Connect 


Ames Research Center 


. XK O — ‘ <S ps 
aerate | ‘. 
one, 


SN 


SF 


‘ 
SJ ww 
Wee 


ort 


Sloan 
Vix 


Total = 25, sparse points away from surface 
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A GRID QUALITY CHECK UTILITIES IN CGT (3) S 
“Bomain Connectivity: Fringe Point Donor Stencil Accuracy 


Histogram of distance between fringe point and vertex obtained by 
donor stencil interpolation (intchk tool in CGT) 


Distance Number of pts % Total 
d<0.0001, 2592370, 89.207 
0.0001 <= d < 0.001, 127886, 4.401 


0.001 <= d <0.01, 128241, 4.413 Set NORFAN carefully 


in OVERFLOW for 


0.01 <=d<0.1, 47312, 1.628 | cue stencil repair 
0.1<=d<1.0, 10167, 0.350 ecclesia ail 
1.0 <=d < 10.0, 49, 1.7E-03 — 


10.0 <=d, 0, 


Donor stencil 


with NORFAN 
25 
point | 
Point obtained using -—- 
interpolation coefficients \—_——— 


Fringe point on 8 corners of donor cell 


a 
-_ ochcenn CID QUALITY CHECK UTILITIES IN CGT (4) w 
Domain Connectivity: Donor Stencil Attributes Compatibility 


Compatibility of cell attributes between fringe point and donor stencil 

- Cell volume ratio histogram table (intchk) and location map (OVERGRID) 
- Bad ratio => gradients cannot be transferred accurately between grids 
Other attributes that could be checked = si 

- Cell aspect ratio, orientation : 


OA 
etremetehatgsn ates 
ober ie 


Volume ratio #Pts. % Tota 


0.5<=C<=1.0 2714268 48.26 


0.2<=C< 0.5 1705036 30.32 
0.1<=C< 0.2 670232 11.92 


0.01<=C< 0.1 525048 9.34 
0.001 <=C< 0.01 9631 0.17 
C< 0.001 21 0.37E-03 


Volume ratio < 0.01 
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x ochcenn CID QUALITY CHECK UTILITIES IN CGT (5) w 
Domain Connectivity: Conversion to Lower Fringe Layers 


- Insufficient grid overlap to support double fringe locally 
- Option to convert from double fringe to single fringe 
=> full 5-point differencing stencil not supported in flow solver 
wer accuracy, robustness) 
pal 
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JAA crip QUALITY CHECK UTILITIES IN CGT (6) wy 
Domain Connectivity: Donor Stencil Quality 


Histogram table (intchk) and location map (OVERGRID) 


SESS 


SSS co 2 1 At 


Stencil Quality Number % Total 


Q=0.0 0 0.0000 
0.0< Q<0.1 0 0.0000 
0.1<=Q<0.2 0 0.0000 
0.2<=Q<0.3 4858 0.1672 
0.3<=Q<0.4 12120 0.4171 
0.4<= Q<0.5 14660 0.5045 
0.5<= Q < 0.6 14054 0.4836 : 
0.6<= Q <0.7 19504 0.6712 
0.7<=Q<0.8 24788 0.8530 
0.8<= Q <0.9 23280 0.8011 
0.9<= Q<1.0 45317 1.5594 
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vy SUMMARY AND CONCLUSIONS (1) Se 
Ames Research Center Workshop Baseline Meshes 


- Grid systems generated and scripted using Chimera Grid Tools 
- Full flap gap geometry (coarse, medium, fine, and extra fine levels) 
- Partially-sealed flap gap (medium only) 


- Workshop guidelines are mostly consistent with current overset grid 
generation best practices 


- Surface grid generation is the most time consuming step 


- Some adjustments needed in developing grid scripts for different levels 
of grid resolution => ideas for further automation development 


- Total development time for all 5 systems ~ 121 man hours 


- Grid script execution time ~ a few minutes (coarse, medium, fine), 
half hr.+ (extra-fine) 


- Preliminary solutions have been computed using OVERFLOW for all 
5 grid systems 
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J SUMMARY AND CONCLUSIONS (2) S 
Grid Quality Checks 


- Effective evaluation using histograms and location maps 


- Wish list 
- Distance to geometry 
- Distance of first volume grid point to wall 
- Cell orthogonality 


- Must-pass grid quality checks 
- Jacobians and self-intersection on surface 


- Mostly-pass grid quality checks 
- Cell volumes 
- Various domain connectivity statistics 
- Need study on how flow solution is affected 
- Accuracy 
- Convergence 
- Robustness / Stability 
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